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1  Introduction 


New  and  improved  methods  of  remote  sensing  have  increased  the  under¬ 
standing  of  Earth’s  origins,  its  resources,  and  those  processes  that  contribute 
to  its  dynamic  (on  a  large  time  scale)  nature.  Scientists  and  engineers  from 
many  disciplines  are  constantly  exploring  new  ways  to  quantify  the  Earth’s 
propenies  for  their  particular  applications.  Among  these  methods  are 
measurements  of  electromagnetic  energy  in  many  different  wavelength 
regimes,  both  passive  and  active.  The  development  of  small  powerful  sources 
and  ultrasensitive  receivers  along  with  improved  data  processing  capabilities 
has  fostered  renewed  interest  in  measurements  of  the  Earth  within  the  micro- 
wave  region  of  the  spectrum  where  wavelengths  in  air  range  from  a  few  milli¬ 
meters  to  several  meters. 

Natural  terrain  surfaces  consist  of  bare  soils,  rocks,  vegetation,  and  water. 
This  study  was  initiated  with  the  hope  of  making  a  meaningful  contribution  to 
the  understanding  of  microwave  interactions  with  namral  terrain.  The  first 
report  in  this  series  attempts  to  summarize  what  is  currently  known  about  the 
measurement  and  modeling  of  the  electrical  properties  of  well-characterized 
soils  (Curtis,  “Microwave  Dielectric  Behavior  of  Soils;  Report  1,  Summary  of 
Related  Research  and  Applications,”  In  Preparation),  while  the  second  report 
contains  a  description  of  a  new  dielectric  property  measurement  capability 
recently  developed  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  in  Vicksburg,  MS  (Curtis,  “Microwave  Dielectric  Behavior  of  Soils; 
Report  2,  A  Unique  Coaxial  Measurement  Apparatus,”  In  Preparation).  It  is 
hoped  that  this  new  measurement  apparatus  will  be  used  in  the  future  to  estab¬ 
lish  a  comprehensive  database  of  soil  properties  covering  a  broad  spectrum  of 
soil  types,  measurement  frequencies,  sample  moisture  contents,  and  sample 
temperatures.  This  final  report  presents  some  initial  measurement  results  and 
discusses  two  modeling  approaches  to  interpreting  those  data.  Most  of  the 
information  contained  in  this  series  of  reports  is  drawn  from  a  doctoral  re¬ 
search  program  recently  completed  by  the  author  (Curtis  1992). 

In  an  effort  to  make  each  report  of  this  series  as  self-contained  as  possible, 
the  rationale  for  pursuing  this  research  program  as  well  as  some  supporting 
information  in  appendices  is  included  within  each  report.  Hopefully,  this  will 
be  useful  for  the  individual  who  might  read  only  one  of  the  reports. 
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Reasons  for  Studying  Electrical  Properties  of  Soils 

Soil  moisture  measurements 

The  dominant  factor  that  controls  the  electrical  behavior  of  soils  is  the 
presence  of  water  (Topp,  Davis,  and  Annan  1980).  Obviously,  one  would 
hope  to  take  advantage  of  this  experimental  fact  to  develop  a  means  of  accu¬ 
rately  and  quickly  measuring  the  moisture  in  soils  without  having  to  collect 
numerous  samples  in  the  field,  weigh  them,  dry  them  for  extended  periods  of 
time,  and  weigh  them  again  to  obtain  either  a  gravimetric  (weight  of  water/ 
weight  of  dry  soil)  moisture  content  or  a  more  useful  volumetric  (volume  of 
water/volume  of  soil  sample)  moisture  content. 

Numerous  attempts  have  been  made  to  develop  a  useful  method  for 
measuring  soil  moisture  content.  All  have  met  with  varying  degrees  of  suc¬ 
cess,  and  none  has  proven  accurate  under  all  conditions.  For  example,  a 
technique  that  essentially  amounts  to  burying  radar  transmit  and  receive 
antennas  in  the  soil  and  relating  the  measurements  of  attenuated  received 
signals  to  moisture  content  (Birchak  et  al.  1974)  is  destructive  to  the  soil 
fabric  (the  way  in  which  soil  particles  are  arranged),  is  very  much  controlled 
by  the  fabric  and  the  size  distribution  of  particles,  and  precludes  the  use  of 
the  same  instrument  in  multiple  locations  (nature  is  not  homogeneous),  as 
well  as  the  ability  to  easily  repair  defective  equipment. 

Another  approach  taken  by  some  researchers  for  making  field  measure¬ 
ments  of  soil  moisture  (that  also  has  numerous  application  in  the  biomedical 
field)  is  that  of  measuring  the  change  in  fringe  capacitance  of  an  open-ended 
coaxial  probe  (Thomas  1966;  Brunfeldt  1987;  Gabriel,  Grant,  and  Young 
1986).  When  pressed  against  a  soil  whose  properties  are  unknown,  the 
resulting  change  in  capacitance  produced  by  the  impedance  mismatch  is 
related  to  electrical  properties  through  calibration  relationships.  Problems 
arising  from  these  measurements  include  the  need  to  have  proper  contact 
between  the  probe  tip  and  the  soil  surface,  the  fact  that  the  volume  of  mate¬ 
rial  associated  with  the  fiinge  capacitance  is  quite  small  (1  cm’  or  less),  and 
that  calibration  conditions  simply  cannot  account  for  all  of  the  dielectric  loss 
mechanisms  that  exist  in  natural  soils.  The  losses  in  moist  soils  can  be 
highly  frequency  dependent  over  a  range  of  several  frequency  decades  on  the 
electromagnetic  spectrum.  A  recent  variation  on  the  open-ended  probe  mea¬ 
surement  scheme  involves  the  use  of  a  waveguide  section  instead  of  a 
coaxial  device  (Parchomchuk,  Wallender,  and  King  1990). 

The  concern  over  small  sample  volumes  can  be  overcome  with  a  redesign 
of  the  open-ended  coaxial  probe  that  replaces  the  solid  outer  conductor  with 
several  pointed  tines  (Campbell  1988)  which  allows  for  the  probe  to  be 
pushed  into  the  surface  of  soft  soils.  The  volume  of  soil  enclosed  by  such 
probes  can  easily  be  tens  of  cubic  centimeters.  The  tined  coaxial  probes  that 
have  been  built  to  date  operate  in  a  frequency  range  that  is  very  much  sub¬ 
ject  to  the  material-dependent  loss  mechanisms  alluded  to  above.  Research 
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continues  that  is  directed  to  better  understanding  these  various  mechanisms 
and  to  fabricate  a  probe  that  operates  in  a  frequency  range  that  is  not  subject 
to  such  material-dependent  anomalies. 


Subterranean  investigations 

There  are  many  types  of  electrical  measurements  in  soils  (and  rocks)  that 
are  used  to  understand  what  lies  beneath  the  Earth’s  surface  (Telford  et  al. 
1984).  Among  these  are  various  schemes  for  measuring  soil  resistivity  as 
well  as  the  attenuation  of  propagating  electromagnetic  waves. 

Resistivity  measurements.  Resistivity  data  in  soils  are  collected  by 
injecting  known  currents  (usually  of  a  frequency  less  than  60  Hz;  alternating 
current  required  to  minimize  effects  of  charge  buildup  on  the  probe)  into  the 
ground  and  measuring  potential  differences  across  pairs  of  nearby  electrodes. 
Assuming  homogeneous  media  and  uniform  resistivity,  it  is  possible  to  calcu¬ 
late  from  the  potential  differences  an  apparent  resistivity  of  the  earth  mate¬ 
rial.  I'hese  measurements  will  be  affected  by  the  presence  of  water  in  the 
soil  or  rock,  by  the  presence  of  mineral  compounds  that  could  go  into  ionic 
solution  with  available  water,  and  by  the  physical  structure  of  the  subsurface 
terrain  itself.  This  being  the  case,  rjsistivity  measurements  are  useful  as  a 
measure  of  subsurface  water  volume,  the  locations  of  mineral  deposits,  and 
subsurface  structure. 

Electromagnetic  wave  propagation.  Another  method  of  making  electri¬ 
cal  subsurface  measurements  involves  the  transmission  of  electromagnetic 
waves  into  the  soil  and  detection  of  energy  that  results  from  waves  reflected 
from  subsurface  anomalies.  Having  the  ability  to  track  wave  propagation  in 
time,  either  by  pulsing  the  source  or  sweeping  over  a  known  frequency  band 
in  some  controlled  manner,  means  that  electromagnetic  wave  propagation 
methods  of  subterranean  investigations  are  particularly  useful  for  locating  the 
depth  of  electrical  anomalies  such  as  the  water  table  in  sandy  soil  (Olhoeft 
1983;  Stewart  1982;  Wright,  Olhoeft,  and  Watts  1984),  buried  pipes  or  wires, 
or  cavities  s.:ch  as  tunnels  or  caves  (Ballard  1983).  Other  applications 
include  the  delineation  of  stratified  media  (Lundien  1972)  and  determination 
of  the  thickness  of  ice  and  frost  layers  (Jakkula,  Ylinen,  and  Tiuri  1980; 
O’Neill  and  Arcone  1991). 

There  are  some  practical  bounds  on  the  utility  of  radio  frequency  systems 
to  conduct  subterranean  investigations  because  of  the  phenomenon  of  "skin 
depth,"  a  measure  of  the  attenuation  of  the  electromagnetic  energy  as  it 
travels  through  the  medium.  For  relatively  low-frequency  sources  (approxi¬ 
mately  200  MHz),  it  can  be  shown  that  low-loss  soils  such  as  dry  sands  can 
possess  a  skin  depth  of  10  to  15  m,  while  high-loss  soils  such  as  wet  silts 
and  clays  may  have  skin  depths  of  only  a  few  centimeters.  Modem  radio 
frequency  receivers  are  extremely  sensitive  devices,  often  having  a  dynamic 
range  of  50  to  100  db.  The  signal  at  skin  depth  represents  about  an  8.7-db 
loss  in  power  or  a  two-way  loss  at  the  receiver  for  reflected  signals  of 
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about  1 7.4  db.  It  is  certainly  not  inconceivable  that  radio  frequency  receiv¬ 
ers  should  be  capable  of  successfully  detecting  reflected  signals  from  subter¬ 
ranean  anomalies  at  depths  of  two  to  three  skin  depths  or  more. 


Remote  sensing  of  environment 

Virtually  all  remote  sensing  of  our  environment  from  airborne  or  space- 
borne  platforms  involves  the  measurement  of  electromagnetic  radiation  from 
the  Earth’s  surface  and/or  atmosphere.  Passive  surveillance  involves  mea¬ 
surements  of  emitted  radiation  and  that  reflected  from  natural  sources  such  as 
the  sun,  the  atmosphere,  and  surrounding  terrain.  Active  remote  sensing 
measurement  systems  include  a  source  for  illuminating  the  target  of  interest. 
Whether  passive  or  active,  visual,  thermal  infrared,  microwave,  or  millimeter 
wave  sensors  are  utilized,  remote  sensing  is  the  collection  and  interpretation 
of  electromagnetic  radiation  and,  as  such,  demands  an  understanding  of  the 
dielectric  properties  of  those  materials  being  observed. 

Environmental  remote  sensing  applications  form  a  list  that  grows  yearly 
as  electronic  components  are  improved  and  data  collection  and  processing 
hardware  and  software  become  faster,  more  reliable,  and  less  costly.  For 
example,  satellites  can  provide  worldwide  surveys  of  land-use  patterns  to 
monitor  the  threat  of  urbanization,  waste  disposal,  and  erosion  of  the  land 
(Colwell  1983).  Similar  systems  (including  those  mounted  in  aircraft)  can 
monitor  the  health  of  vegetation  to  keep  abreast  of  such  things  as  loss  of 
forest  and  the  potential  for  food  shortages.  Sea  traffic  in  the  far  northern  and 
southern  shipping  lanes  can  be  made  safer  through  the  use  of  airborne  and 
spacebome  sensors  to  detect  ice  hazards. 

One  of  the  more  obvious  applications  for  microwave  remote  sensing 
devices  is  that  of  conducting  surface  moisture  surveys  to  help  predict 
groundwater  availability  and  the  potential  for  flooding.  Attempts  have  been 
made  to  relate  soil  moisture  to  both  laboratory  reflectance  data  (Lundien 
1966)  and  to  airborne  sensor  backscatter  measurements  (Ulaby,  Cihlar,  and 
Moore  1974;  John  1992).  Careful  airborne  sensor  measurements  might  pro¬ 
vide  a  first  approximation  to  the  complex  dielectric  constant  of  the  soil  near 
the  surface. 

Because  of  the  difference  in  dielectric  behavior  of  liquid  water  and  vari¬ 
ous  forms  of  ice,  it  may  be  possible  to  use  airborne  sensors  to  detect  freezing 
and  thawing  (Wegmuller  1990)  in  remote  locations  that  could  be  used  to 
predict  spring  runoff  conditions  and  all  that  that  encompasses  for  agricultural 
applications,  the  effects  on  the  fishing  industry,  and  the  anticipation  of  flood¬ 
ing  in  built-up  areas.  Military  analysts  are  concerned  about  soil  moisture 
conditions  because  of  its  impact  on  trafficability,  the  ability  of  vehicles  to 
move  effectively  over  natural  terrain. 

Another  remote  sensing  application  is  the  mapping  of  exposed  soils  and 
rocks  in  remote  areas  of  the  world  from  high-flying  aircraft  or  satellites. 
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which  might  prove  useful  for  geomorphoiogical  studies  (Swanson  et  al. 

1988)  or  even  mineral  exploration.  A  less  obvious,  but  recent,  application  of 
microwave  remote  sensing  in  soils  dealt  with  archeological  surveys  in  desert 
areas  (Berlin  et  al.  1986;  McCauley  et  al.  1986) 


Others 

While  the  above  paragraphs  emphasize  some  of  the  most  obvious  and 
useful  applications  of  a  better  understanding  of  soil  electrical  properties, 
others  have  been  noted  in  the  literature.  For  example,  some  researchers  have 
attempted  to  relate  electrical  property  measurements  to  the  physical  proper¬ 
ties  of  soils  (Campbell  and  Ulrichs  1969;  Hayre  1970;  Arulanandan  and 
Smith  1973;  Madden  1974).  Of  course,  nothing  has  been  said  about  the 
military’s  need  to  better  understand  the  microwave  response  of  soils  that 
form  the  backgrounds  to  military  targets;  i.e.,  when  and  why  does  clutter 
become  a  source  of  target-like  signatures?  Outside  of  the  topic  of  soils, 
studies  of  the  microwave  response  of  foodstuffs  has  direct  application  to 
quality  control  concerns  in  the  food  industry  (Nelson  1973,  1983). 

Another  new  application  of  this  technology  that  is  closely  related  to  the 
discussion  on  soil  moisture  is  that  of  detecting  liquid  ground  contaminants, 
either  near  the  surface  or  at  arbitrary  depths  using  a  specially  fabricated 
probe.  If,  as  will  be  argued  later,  polarizable  liquids  can  be  characterized  by 
a  unique  frequency  of  peak  losses  because  of  the  dielectric  relaxation  phe¬ 
nomenon,  then  a  probe  could  be  designed  to  measure  losses  over  a  frequency 
span  that  is  broad  enough  to  detect  a  peak  loss  frequency  and,  coupled  with 
the  results  of  a  thorough  experimental  program,  to  identify  the  particular 
contaminant. 


Problem  Statement 

The  overall  objective  of  this  research  is  to  measure  and  model  the  dielec¬ 
tric  response  of  moist  soils  in  ways  that  will  test  the  current  understanding  of 
loss  mechanisms  over  a  broad  range  of  moisture  levels,  frequencies,  and 
material  temperatures  and  perhaps  suggest  new  ways  of  thinking  about  losses 
as  a  function  of  these  variables. 

The  first  step  in  achieving  this  objective  has  already  been  accomplished 
through  the  development  of  an  apparatus  that  allows  for  the  measurement  of 
attenuation  in  moist  soils  over  a  broad  spectrum  of  frequencies  at  controlled 
temperatures  and  calculable  moisture  contents  (Curtis,  “Microwave  Dielectric 
Behavior  of  Soils;  Report  2,  A  Unique  Coaxial  Measurement  Apparatus,”  In 
Preparation).  The  remainder  of  this  report  describes  measurements  made  on 
a  number  of  different  soils  and  two  approaches  to  modeling  the  measured 
response. 
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Sample  Preparation  and  Measurements 

Soil  samples  were  prepared  from  one  of  four  distinctly  different  soils. 
These  consisted  of  a  poorly  graded  sand,  a  well-graded  sand,  a  poorly  graded 
clean  silt,  and  almost  pure  kaolinite,  a  nonswelling  clay  mineral.  The  soils 
and  their  associated  properties  were  obtained  from  the  Geotechnical  Labora¬ 
tory  at  the  WES.  Each  soil  type  batch  was  made  as  homogeneous  as  possi¬ 
ble  by  drying,  pulverizing,  and  mixing  thoroughly.  Samples  were  then  taken 
to  determine  the  usual  engineering  properties  of  grain  size  distribution  and 
plasticity  indices.  The  results  of  those  measurements  are  contained  in 
Appendix  A,  along  with  the  results  of  specific  surface  measurements  on  each 
sample  that  were  conducted  by  a  commercial  laboratory. 

The  normal  measurement  procedure,  as  previously  described  in  Report  2 
of  this  series,  was  the  following: 

a.  The  dry,  empty  sample  holder  was  carefully  weighed. 

b.  A  nearly  saturated  sample  of  soil  was  prepared  using  distilled,  deion¬ 
ized  water  (a  reasonable  amount  of  time  was  allowed  for  the  silt  and 
clay  to  reach  some  sort  of  equilibrium)  and  was  placed  in  the  sample 
holder.  The  holder  was  tapped  against  a  flat,  hard  surface  to 
uniformly  settle  the  sample. 

c.  The  sample  and  holder  were  carefully  weighed,  then  sealed  and  con¬ 
nected  to  the  measurement  system. 

d.  The  temperature  of  the  bath  was  set  at  -10  “C.  For  each  new  temper¬ 
ature  setting,  the  sample  was  given  10  min  to  reach  equilibrium  once 
the  bath  had  reached  the  desired  setting.  The  10-min  wait  was  found 
experimentally  to  be  more  than  adequate  for  sample  equilibrium  to  be 
achieved. 

e.  The  fi^uency  range  for  the  dual  directional  coupler  being  used  was 
stepped  through,  and  calculations  of  e'  and  e"  were  made. 
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/  The  temperature  of  the  sample  was  changed  and  measurements 
repeated  in  the  following  sequence:  -10,  -5.  -2,  0,  2,  5,  10,  20,  30, 
and  40  °C. 

g.  The  holder  was  unsealed  and  its  cover  loosened.  The  sample  and 
holder  were  periodically  weighed  to  determine  approximate  values  of 
moisture  content.  When  the  sample  was  thought  to  be  at  about  the 
desired  volumetric  moisture  content,  the  holder  was  resealed  and 
reconnected;  the  sample  temperature  was  taken  back  to  -10  °C.  and 
the  whole  measurement  sequence  repeated. 

h.  After  all  measurements  were  conducted,  the  sample  was  carefully 
removed,  dried,  and  weighed  to  detet<  i:ne  its  original  dry'  density  and 
to  facilitate  calculations  of  volumetric  moisture. 

As  explained  in  Report  2,  measurements  had  to  be  conducted  over  two 
different  frequency  ranges  because  of  the  unavailability  of  a  single  dual 
directional  coupler  that  could  cover  the  entire  range  desired. 


Summary  of  Data  Collected 

Figure  1  serves  as  an  indication  of  both  data  quantity  and  quality.  On 
this  figure  are  symbols  that  represent  data  collected  for  the  four  primary  soil 
materials  at  one  temperature,  20  °C.  For  most  of  these  data  points,  data 
were  collected  at  nine  other  sample  temperatures.  Although  they  are  not 
shown  on  this  figure,  data  were  also  collected  for  the  empty  holder  and  for 
water,  alcohols,  and  a  swelling  clay  mineral  called  hectorite.  Results  for 
hectorite  are  not  reported  because  physical  property  information  is  not  cur¬ 
rently  available  for  that  material.  The  total  number  of  data  sweeps  (a  sweep 
of  frequencies)  collected  during  these  experiments  easily  exceeds  five 
hundred. 

The  reference  above  to  data  quality  is  in  terms  of  the  repeatability  of 
sample  dry  density  as  a  result  of  the  crude  sample  preparation  technique 
described  above.  The  strong  clustering  of  the  data  points  for  each  material 
type  indicates  how  closely  sample  dry  density  matched  for  all  samples.  Only 
one  string  of  data  is  significantly  shifted  from  the  others,  and  that  is  the  data 
for  kaolinite.  The  four  data  points  that  lie  to  the  right  of  the  others  represent 
a  sample  with  a  dry  density  approximately  20  percent  greater  than  the  other 
clay  mineral  samples.  The  starting  point  for  this  sample  was  a  batch  of  kaol¬ 
inite  that  had  been  allowed  to  dry  somewhat  before  the  sample  was  inserted 
into  the  sample  holder. 
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Figure  1 .  Data  collected  at  20  ”C 

Dispersion  in  Soils  as  a  Function  of  Moisture 
Content 


If  the  electrical  behavior  of  moist  soils  is  controlled  by  the  amount  and 
purity  of  water  and  how  it  fills  the  interstices,  then  one  should  be  able  to 
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observe  high-frequency  dielectric  relaxation  and  low-frequency  losses 
because  of  conductivity  or  Maxwell- Wagner  effects  or  bound  water  (Curtis. 
“Microwave  Dielectric  Behavior  of  Soils;  Report  I,  Summary  of  Related 
Research  and  Applications,”  In  Preparation),  and  these  phenomena  should  be 
more  pronounced  as  moisture  content  increases.  One  should  also  observe 
different  responses  in  different  soils  at  the  same  moisture  content  because  of 
such  factors  as  specific  surface  and  the  ways  in  which  water  molecules  are 
attracted  to  the  solid  surfaces. 

Figures  2-6  summarize  dispersion  measurements  for  the  four  soil  types. 
Each  figure  contains  data  for  sample  temperatures  of  either  -10  °C  or 
+20  °C.  These  two  temperatures  were  chosen  simply  to  compare  soils  with 
frozen  water  and  soils  with  liquid  water.  The  vertical  dashed  line  at  2  GHz 
represents  the  break  point  between  measurements  made  with  the  low- 
frequency  apparatus  configuration  and  those  made  with  the  high-frequency 
setup.  Three  figures  contain  silt  data,  the  third  one  representing  results  for  a 
slightly  different  experimental  procedure.  In  the  latter  case,  both  frequency 
range  sweeps  were  conducted  before  the  sample  was  incrementally  dried.  It 
is  very  satisfying  that  the  results  for  the  two  setups  overlap  almost  exactly. 

Conclusions  that  may  be  drawn  from  these  data  include  the  following; 

a.  Frozen  soils  do  not  exhibit  high-frequency  dielectric  relaxation,  but  do 
reveal  some  low-frequency  losses. 

b.  Nonfrozen  soils  do  demonstrate  high-frequency  relaxation,  which  does 
increase  with  increasing  moisture  content.  These  results  are  quite 
consistent  with  those  published  by  Hallikainen  et  al.  (1985). 

c.  Low-frequency  losses  in  nonfrozen  soils  are  also  proportional  to  mois¬ 
ture  content.  Values  are  consistent  with  those  of  Campbell  (1988). 
Minimum  losses  in  nonfrozen  soils  occur  in  the  1-  to  2-GHz  range, 
which  agrees  with  the  observations  of  Hallikainen  et  al.  (1985). 

d.  Comparison  of  sand,  silt,  and  clay  high-frequency  losses  at  compara¬ 
ble  moisture  contents  provides  no  indication  of  differences  in  free 
water  content  in  spite  of  the  great  variation  in  specific  surfaces. 

e.  If  the  relative  magnitudes  and  frequency  ranges  of  different  loss 
mechanisms  published  by  Hasted  (1973)  are  correct,  then  bound  water 
relaxation  is  a  good  candidate  for  some  of  the  low-frequency  losses 
seen  in  these  data.  One  cannot,  however,  discount  ionic  conductivity 
or  Maxwell- Wagner  losses,  or  both,  as  contributing  loss  mechanisms. 

/  The  electrical  response  of  frozen  soils  appears  to  be  relatively  insensi¬ 
tive  to  moisture  content  within  the  frequency  range  measured  in  this 
report,  but  does  show  some  dependence  on  soil  type  at  low 
frequencies. 
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Figure  5.  Dispersion  in  silt  (second  set) 
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From  the  results  of  their  excellent  measurements  on  real  soils  ranging 
from  sandy  loams  to  heavy  clays.  Topp,  Davis,  and  Annan  (1980)  reported 
that,  in  the  frequency  range  of  20  MHz  to  1  GHz,  the  apparent  permittivin 
“was  strongly  dependent  on  (moisture  content)  and  only  weakly  on  soil  type, 
density',  temperature,  and  frequency.”  While  the  measurements  reported  in 
this  study  have  not  dealt  with  density  effects,  and  they  do  support  the  con¬ 
tention  that  temperature  effects  are  minimal  for  unfrozen  soils,  they  clearly 
demonstrate  that  soil  type  and  frequency  are  important  variables. 


Single  Frequency  Observations 

Temperature  effects 

Data  collected  at  -10,  -5,  -2,  0,  2,  5,  10,  20,  30,  and  40  °C  are  quite  ade¬ 
quate  for  visualizing  the  temperature-dependent  behavior  of  the  complex 
dielectric  constant  at  selected  frequencies.  These  results  are  shown  on  Fig¬ 
ures  7-10  and  lend  themselves  to  several  observations; 

a.  Liquid  water  appears  at  temperatures  less  than  0  °C,  probably 
somewhere  between  -2  and  -5  °C.  This  is  consistent  with  results 
reported  by  Delaney  and  Arcone  (1982). 

b.  The  real  part  of  the  dielectric  constant  for  nonfrozen  soils  is  approxi¬ 
mately  independent  of  temperature  for  ail  soil  types. 

c.  The  imaginary  part  of  the  dielectric  constant  for  nonfrozen  soils  is  not 
independent  of  temperature,  showing  a  rise  with  temperature  at  low 
frequencies  and  a  drop  with  temperature  at  high  frequencies. 

Because  the  low-frequency  data  show  a  rise  in  the  loss  term  with  increas¬ 
ing  temperature,  an  estimate  of  the  activation  energy  involved  in  the  loss 
mechanism  can  be  made  following  the  method  of  Hoekstra  and  Doyle  (1971) 
and  Campbell  (1988),  which  was  discussed  earlier  in  Report  1.  Using  their 
formula 


ln(e")  =  - 


A 


1 


Kelvin 


(1) 


where 


e"  =  imaginary  part  of  the  complex  relative  dielectric  constant 

A  =  activation  energy,  or  the  energy  needed  to  overcome  some 
equilibrium  state 
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R  =  universal  gas  constant 

^Kelvin  ~  absolute  temperature 

and  the  temperature  effects  data  at  100  MHz  yields  an  average  (over  mois¬ 
ture  contents)  activation  energy  of  3.1  Kcal/mole  for  silt  and  3.3  Kcal/mole 
for  kaolinite.  These  numbers  are  within  the  range  of  values  attributed  to 
hydrogen  bond  breakage. 

Moisture  effects 

Figures  11-14  summarize  data  at  selected  frequencies  and  at  +20  °C  as  a 
function  of  volumetric  moisture  content  for  each  of  the  materials  tested. 

What  can  be  said  immediately  from  an  examination  of  the  data  on  moist 
soils  presented  in  this  format  is  the  following; 

a.  The  real  part  of  the  dielectric  constant  for  all  nonfrozen  soils  is  best 
described  by  some  nonlinear  relationship  with  volumetric  moisture. 
Bilinear  fits  such  as  those  sketched  in  on  each  figure  would  break  at 
critical  values  of  moisture  content  in  the  range  of  30  to  35  percent  for 
the  silt  and  clay  minerals  and  in  the  10-  to  15-percent  range  for  the 
sands.  Wang  and  Schmugge  (1980)  reported  critical  moisture  con¬ 
tents  in  the  20-  to  30-percent  range  for  the  real  soils  that  they  studied. 

b.  The  capacitive  nature  of  nonfrozen  soils  as  reflected  in  the  real  part  of 
the  dielectric  constant  varies  with  soil  type. 
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Figure  14.  Moisture  effects  for  nonfrozen  soils  at  8  GHz,  20  ®C 
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3  Equivalent  Circuit  Modeling 


The  following  paragraphs  describe  an  effort  to  analytically  model  the  mea¬ 
sured  soil  response  through  an  adaptation  of  the  classical  technique  of  repre¬ 
senting  the  electrical  response  of  an  heterogeneous  mixture  by  an  equivalent 
circuit.  The  models  developed  here  incorporate  the  concept  of  series  and 
parallel  behavior  discussed  in  Report  1  of  this  series,  as  well  as  a  model  for 
water  behavior  that  accounts  for  dielectric  relaxation.  The  new  features  of 
this  model  are  the  broad  frequency  band  of  simulations  through  a  modification 
of  the  equivalent  circuit  for  water  and  the  attempt  to  fix  the  values  of  as  many 
circuit  elements  as  possible  while  giving  physical  interpretations  to  the  remain¬ 
ing  parameters. 


Equivalent  Circuit  Representation 

If  one  accepts  that  the  electrical  response  of  soil-water  mixtures  is 
bounded  on  the  lower  end  of  moisture  contents  by  a  "series-like"  behavior 
and  on  the  upper  end  of  moisture  contents  by  a  "parallel-like"  behavior 
(Campbell  1988),  then  the  most  sensible  simulation  should  include  both  ele¬ 
ments.  With  the  requirement  for  both  series  and  parallel  response,  one  sound 
approach  to  modeling  the  response  of  moist  soil  by  the  method  of  equivalent 
circuits  is  the  three-paUi  model  first  suggested  by  Sachs  and  Spiegler  (1964) 
and  adapted  by  Smith  (1971).  However,  the  specific  representation  used  by 
Smith  (see  Figure  36  in  Report  1)  is  not  adequate  to  simulate  both  low- 
frequency  loss  mechanisms  and  losses  at  higher  frequencies  due  to  dielectric 
relaxation.  A  variation  of  the  three-path  model  was  adopted  for  this  study 
using  the  following  reasoning. 

Data  from  the  last  chapter  and  from  other  sources  show  clearly  that  the 
response  of  moist  soils  is  controlled  by  the  presence  of  water  in  the  samples. 
Most  of  the  data  show  a  high-frequency  anomalous  loss  highly  correlated 
with  moisture  content  and  a  conductivity-like  low-fi^uency  loss  mechanism. 
In  fact,  the  data  do  not  say  whether  the  low-fi^uency  losses  are  due  to  a 
conductivity  effect  in  free  water  or  a  Maxwell-Wagner  mechanism  caused  by 
the  presence  of  the  soil  particles  as  suggested  by  many.  Lower  frequency 
measurements  coupled  with  a  much  more  thorough  understanding  or  control 
of  the  soil  chemistry  would  help  answer  this  last  concern.  Nevertheless,  the 
best  equivalent  circuit  for  the  liquid  paths  has  to  be  one  that  includes  a 
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high-frequency  permittivity,  a  direct  current  conductivity,  and  relaxation 
losses  in  the  10-  to  20-GHz  range  of  frequency  (see  Figure  34,  Report  1). 

As  for  the  contribution  by  the  soil,  a  good  representation  for  its  electrical 
response  is  felt  to  be  one  that  accounts  for  a  constant  wave  speed  at  all  fre¬ 
quencies  while  including  a  low-frequency  loss  term.  Therefore,  a  parallel 
circuit  (see  Figure  30  (b).  Report  1)  was  chosen  as  the  most  appropriate 
representation  for  the  soil  component  of  the  mixture.  One  implication  of  this 
selection  is  that  one  cannot  model  anomalous  losses  within  the  soil  that  result 
from  the  interaction  of  the  sc:'  and  water  particles.  However,  in  the  range  of 
frequencies  covered  by  these  experiments,  it  is  probably  not  possible  to 
detect  the  difference  in  losses  because  of  one  mechanism  or  the  other. 

Combining  these  two  sets  of  elements,  the  final  three-path  equivalent 
circuit  used  to  simulate  the  response  of  the  soils  tested  in  this  study  is  that 
shown  in  Figure  15,  where  the  parameters  “a,”  “b,”  “c,”  and  “d”  are  weight¬ 
ing  constants  whose  physical  significance  will  be  discussed  later.  A  thought 
process  for  electing  the  model  parameters  is  discussed  later,  and  the  results 
of  a  number  of  simulations  will  be  shown;  but  first,  a  brief  presentation  of 
the  governing  equation  for  this  circuit  is  given. 


Figure  15.  Three-path  equivalent  circuit  used  for  this  study 

The  macroscopic  complex  dielectric  constant  of  the  moist  soil  material  is 
modeled  as  the  equivalent  capacitance  of  the  circuit  shown  in  Figure  15. 
Using  the  w  subscript  for  water  and  the  s  subscript  for  soil,  the  macroscopic 
dielectric  properties  can  be  written 
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e  =  e'  +  it"  =  oE 


+  +  C£ 


(2) 


where 


(3) 


In  terms  of  the  simple  circuit  capacitances  and  resistances. 
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Parameter  Selection  and  Model  Execution 

A  reasonable  set  of  circuit  parameters  begins  with  the  same  values  chosen 
for  water  that  are  shown  on  Figure  34  of  Report  1. 

C.,  =  4 


C.,  =  76 


R^i  =  5.888x10'*  sec/cm 
-  l-25xl0''^  sec/cm 

They  allow  for  a  dielectric  relaxation  peak  at  about  16.75  GHz  and  for  finite 
values  of  conductivity  losses  at  frequencies  less  than  1  MHz.  Although  the 
data  measured  in  this  study  did  not  go  as  low  as  1  MHz,  it  was  discovered 
by  trial  and  error  that  the  low-frequency  conductivity  term  helped  improve 
the  model’s  fit  to  data  through  its  interaction  with  the  soil  element  response. 
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As  for  the  soil  parameters,  a  static  permittivitv  C,  was  taken  to  be  about 
2.  That  only  leaves  the  soil  resistivity  to  be  considered  It  so  happens 
that  this  model  allows  soil  resistivity  to  be  a  variable;  however,  it  is  rather 
easy  to  estimate  its  range  in  the  following  way.  The  model  has  associated 
the  low-  frequency  losses  with  the  soil  element.  They,  in  turn,  are  l/(a)R‘). 

A  typical  value  of  the  loss  term  is  10  at  2  GHz  for  wet  silt.  This  leads  to  an 
estimate  of  the  resistivity  as  about  8  x  10  '*  sec/cm. 

To  obtain  values  of  the  weighting  parameters,  “a,’'  “b.’"  “c.”  and  “d."  as 
well  as  the  value  of  /?,,  a  simple  iterative  code  was  developed  that  minimized 
the  difference  between  circuit  predictions  and  real  data.  Smith  (1971)  used  a 
similar  approach  to  calculate  the  parameters  associated  with  his  three-path 
model  and  hinted  that  it  was  a  rather  sophisticated  calculation  process.  In 
truth,  it  is  not  necessary  to  do  anything  more  complicated  than  to  increment 
each  parameter  of  the  equivalent  circuit  (element  values  and  weighting 
parameters)  and  to  calculate  the  difference  between  all  data  points  and  pre¬ 
dicted  values  for  each  parameter  change  taken  one  at  a  time.  That  parameter 
increment  that  caused  the  greatest  reduction  in  difference  was  applied  at  each 
step.  Parameter  values  were  both  increased  and  decreased  for  each  iteration, 
and  the  magnitude  of  these  changes  were  decreased  as  the  model  simulation 
converged. 

Smith  (1971)  also  chose  to  allow  all  of  the  parameters  in  his  three-path 
model  to  vary  except  for  the  permittivity  of  water.  In  this  model,  most  of 
the  electrical  parameters  were  fixed  because  they  have  physical  meaning. 

The  permittivities  selected  for  both  the  water  element  and  the  soil  element 
reflect  observed  values  that  are  associated  with  wave  velocity  in  each  media. 
The  value  of  resistivity  in  the  water  element  that  results  in  the  high 
frequency  anomalous  behavior  was  selected  to  properly  locate  the  peak 
losses  due  to  dielectric  relaxation.  In  the  absence  of  data  at  higher  frequen¬ 
cies,  the  frequency  of  peak  loss  was  taken  to  be  constant  for  all  moisture 
contents.  The  soil  resistivity  was  allowed  to  be  variable  to  help  optimize  the 
fits  to  low-frequency  loss  data. 


Simulation  Results 


Figures  16-18  contain  a  comparison  of  model  fits  to  real  data  for  tan 
sand,  silt,  and  kaolinite  samples  at  three  different  volumetric  moisture  condi¬ 
tions.  Besides  demonstrating  the  ability  to  simulate  the  real  response  of  the 
soil  samples  quite  well,  it  is  gratifying  to  note  sensible  trends  in  the  weight¬ 
ing  parameters.  Obviously,  as  moisture  levels  rise,  the  contribution  to  sam¬ 
ple  response  from  the  parallel  (or  long-range  connectivity)  water  element 
should  increase  as  reflected  in  an  increase  in  the  parameter  “c.”  Because 
water  dominates  the  electrical  response,  the  soil-related  parameter  “b”  should 
drop  as  moisture  increases.  Similarly,  if  the  series  elements  can  account  for 
Maxwell- Wagner  like  behavior  in  the  soil-water-air  mixture,  then  it  is  also 
sensible  for  the  losses  contributed  by  the  soil  particles  in  the  presence  of  an 
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(a)  7  percent  moisture 


Figure  16.  Equivalent  circuit  model-data  comparisons  for  tan  sand  (Sheet  1  of  3) 
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(c)  27  percent  moisture 


Figure  16.  (Sheet  3  of  3) 
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(b)  25  percent  inoisture 


Figure  17.  (Sheet  2  of  3) 
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(a)  3.5  percent  moisture 


Figure  18.  Equivalent  circuit  model-data  comparisons  for  clay  (Sheet  1  of  3) 
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(b)  42  percent  moisture 


2 

4 

76 

7.253B305G4BE-t0 

5.B08B88E-8 

1.25E-13 

.7V2059299585 

. 157215804939 

.0738222377229 


Chapter  3  Equivalent  Circuit  Modeling 


41 


e'  *■ 


0000  laboratory  data 
-  model  Tit 


1  i« 

FREQUENCY  tCHl) 


5.  19262 I02B9SE-1 1 

5.8880B8E-8 

1.25E-13 

.52899037522 

.451648706606 

.  120798813955 


(c)  64  percent  moisture 


Figure  18.  (Sheet  3  of  3) 
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increasing  amount  of  water  to  increase  as  well.  This  would  be  reflected  in 
an  increase  in  the  parameter  “d”  that,  in  general,  is  observed. 

While  the  arrangement  of  elements  in  the  three-path  model  is  physically 
reasonable  and  the  simulations  are  very  accurate,  a  number  of  questions  con¬ 
cerning  its  application  remain  unanswered.  For  example,  when  does  one 
hold  circuit  element  parameter  values  constant,  and,  if  not,  how  should  those 
parameters  be  allowed  to  vary?  It  might  seem  preferable,  at  first  thought,  to 
hold  all  circuit  element  values  fixed  and  allow  only  the  weighting  constants 
to  vary.  In  fact,  this  approach  was  tried,  but  with  only  limited  success.  If 
only  the  weighting  coefficients  could  vary,  then  they  would  have  to  account 
for  all  of  the  dramatic  changes  that  take  place  chemically  within  the  mixture, 
and  they  do  not  appear  to  be  able  to  do  so.  In  other  words,  by  allowing  only 
weighting  coefficients  to  change  implies  that  the  properties  of  the  basic 
constituents  do  not  change  as  the  mixture  changes.  Even  if  the  selected 
parameters  for  dry  soil  particles  and  unbound  liquid  water  are  quite  accurate, 
weighting  parameters,  by  themselves,  cannot  be  expected  to  account  for  the 
loss  mechanisms  because  of  bound  water  or  because  of  Maxwell-Wagner 
type  effects,  or  enhanced  conductivity  because  of  salts  going  into  solution,  or 
enhanced  capacitive  response  set  up  by  the  parallel  plate-like  structure  of  the 
clay  particles  being  filled  by  water  in  the  interlayer  spaces  or  simply  between 
particles  in  nonswelling  clays. 

If,  then,  some  of  the  circuit  parameters  are  allowed  to  vary,  which  ones 
should?  Why  let  the  resistivity  of  the  soil  element  vary  and  not  the  water? 

A  number  of  simulations  were  tried  with  /?,  fixed  and  free.  These  results 
were  unsatisfactory  for  the  mixtures  that  showed  increasing  £'  values  with 
decreasing  frequency,  with  the  simulation  values  of  e'  becoming  constant 
with  decreasing  frequency. 

In  summary,  the  three-path  equivalent  circuit  is  a  promising  tool  for 
exploring  the  complex  dielectric  behavior  of  heterogeneous  mixtures.  Simu¬ 
lations  done  in  this  study  very  accurately  reproduce  measured  responses. 

The  model  accounts  for  both  series  and  parallel  electrical  responses  and  pos¬ 
sesses  some  parameters  that  can  be  fixed  to  values  measured  in  other  experi¬ 
ments.  There  are  still,  however,  unresolved  issues  such  as  what  physical  or 
chemical  properties  can  be  associated  with  the  parameters  that  do  vary. 

These  issues  should  provide  a  rich  opportunity  for  further  research. 
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4  Fractal  Geometry  Model  and 
Critical  Water  Content 


Another  analysis  technique  used  in  this  study  involves  the  application  of 
fractal  geometry  (Appendix  B)  to  model  the  structure  of  the  porous  media  and 
to  associate  a  change  in  the  fractal  dimension  with  the  onset  of  long-range 
connectivity  .  While  this  model  has  been  used  in  basic  research  on  the 
properties  of  coal,  this  is  its  first  application  to  soils.  The  proposed  hypothe¬ 
sis  is  that  the  soil  moisture  content,  or  cumulative  pore  volume,  at  the  limit  of 
fractal  behavior  is  equivalent  to  the  critical  moisture  that  defines  a  transition 
from  series  to  parallel  electrical  response.  This  hypothesis  was  tested  and 
yielded  fractal  dimensions  for  the  soils  tested  within  this  study  that  are  sup¬ 
ported  by  other  references. 

Campbell  (1988)  suggested  that  water  may  fill  the  pore  spaces  of  a  soil 
sample  in  a  fractal  manner.  Pursuit  of  this  concept,  particularly  as  it  relates 
to  the  distribution  of  pore  sizes  in  a  soil  sample,  has  led  to  the  discovery  of 
related  research  in  pore-size  distributions  (Sridharan,  Altschaeffl,  and 
Diamond  1971;  Arya  and  Paris  1981)  and  in  the  modeling  of  porous  media  as 
fractal  geometries  (Friesen  and  Mikula  1987)  that  lends  itself  to  a  relatively 
simple  fractal  interpretation  of  pore-size  distribution  data.  The  difference 
between  the  following  study  and  Campbell’s  model  is  that  rather  than  assum¬ 
ing  pore-filling  water  clusters  have  a  fractal  geometiy,  this  new  approach 
models  the  soil  structure,  itself,  using  fractal  concepts  and  makes  certain 
assumptions  to  relate  volumetric  soil  moisture  content  to  the  fractal  model 
behavior. 


Fractal  Model  of  Pore-Size  Distribution 


The  porous  soil-water-air  structure  can  be  described  as  a  sponge-like 
structure  similar  to  that  shown  in  Figure  B3.  Much  of  what  follows  has  been 
adopted  from  Friesen  and  Mikula  (1987),  where  additional  references  are 
given  to  some  earlier  work.  Their  application  was  toward  a  better  way  of 
quantifying  the  porous  structure  of  coal. 
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Consider,  first,  a  unit  solid  cube  of  soil  minerals.  Divide  the  cube  equally 
into  subcubes  of  size  d  =  I/m.  Next,  remove  the  subcubes  necessary  to  give 
the  basic  structure  generator  shape.  (In  the  case  of  the  Menger  sponge 
shown  in  Figure  B3,  m  =  3,  and  the  number  of  subcubes  removed  is  7.) 

That  leaves  a  number  of  remaining  cubes  =  A'^.  But  by  the  definition  of 
fractals  given  in  Equation  B2, 


K  =  (d) 


~D 


(6) 


where  D  is  the  fractal  "dimension"  of  the  porous  media.  Now'  repeat  the 
entire  process  on  each  subcube.  After  this  second  step,  the  new  subcubes 
have  a  dimension  of  d‘  and  the  number  of  remaining  subcubes  is  N^.  Con¬ 
tinuing  to  the  hh  subdivision,  there  is  a  particle  size  of  and  the  number  of 
remaining  subcubes  is  N^.  In  other  words,  the  volume  of  remaining  solid 
material  at  step  k  is 


Noting  that  at  each  step,  the  particle  size  is  identical  to  the  smallest  pore 
size,  and  renaming  the  particle  size  (d^  as  l^, 


V,  .  (;,)>■”  (8) 

and  the  pore  volume  is  just 

- (i.r 


What  one  now  has  is  a  fractal  model  for  a  porous  medium  that  relates 
pore  volume  (which  has  to  be  related  to  moisture  content)  and  pore  size. 
The  fractal  dimension  of  the  porous  media  can  be  calculated  by  plotting  the 
slope  of  the  distribution  versus  pore  size  on  a  log-log  plot: 


log 


dy 

port 

~dr 


=  constant 


(2-Z?)log(/) 


(10) 
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Fractal  Model  and  Pressure  Plate  Data 


Unfortunately,  simple  pore-size  distribution  data  do  not  exist.  What  is 
available  is  a  limited  amount  of  porosimetry  and/or  pressure  plate  data  that 
can  be  interpreted  as  pore  distribution  data  through  the  capillary  relationship 
(Bear  and  Verruijt  1987) 


P  = 


cos(0) 

D 


(11) 


where 

p  =  equilibrium  pressure 
6  =  surface  tension  of  nonwetting  fluid 
6  =  contact  angle  between  fluid  meniscus  and  pore  surface 
D  =  capillary  diameter 

Substituting  the  capillary  relationship  into  the  fractal  model  for  pore 
distribution, 

'  =  1  -  K'p^-^  (12) 

where  K  and  K"  represent  a  lumping  of  constants.  Then 

=  -K'{D  -  'i)p°-*  (13) 

dp 

from  which 


log 


dp 


=  constant  -  {D  -  4)log(p) 


(14) 


Therefore,  the  slope  of  a  log-log  plot  of  the  derivative  of  a  porosimetry  or 
pressure  plate  curve  is  proportional  to  D  -  4. 

As  a  test  of  this  relationship,  real  pressure  plate  data  from  Arya  and  Paris 
(1981)  for  five  different  soil  samples  were  digitized,  their  slopes  estimated 
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by  difference  methods,  and  the  results  plotted  on  Figure  19.  In  general,  these 
curves  reveal  a  bilinear  response  with  the  break  point  occurring  at  approxi¬ 
mately  the  field  capacity  of  0.333  bars.  The  field  capacity  is  associated  with 
the  amount  of  water  that  would  drain  from  a  soil  sample  due  to  gravitational 
forces  only.  At  pressures  higher  than  field  capacity  but  less  than  15  bars, 
which  soil  scientists  refer  to  as  the  wilting  point  (Wang  and  Schmugge  1980), 
a  straight  line  fit  to  the  data  on  Figure  19  would  result  in  a  fractal  dimension 
for  the  soil  structure  of  approximately  2.67.  The  break  in  the  slope  of  Lhese 
data  roughly  corresponds  to  volumetric  moisture  contents  of  25  to  40  percent. 


An  interesting  question  is  whether  or  not  one  might  be  able  to  predict  these 
observed  moisture  contents  from  ancillary  data.  For  example,  if  the  fractal 
behavior  of  the  soil  is,  in  fact,  associated  with  water  being  held  by  capillary 
forces,  then  the  field  capacity  can  be  identified  with  a  certain  pore  diameter 
through  the  capillary  equation.  Taking  conditions  for  pure  water  and  clean 
glass  (contact  angle  equal  zero  degrees  and  surface  tension  equal  to 
72  dynes/cm  at  25  °C)  as  an  approximation  for  the  elements  in  moist  soil,  one 
can  calculate  from  the  capillary  equation  that  the  range  of  pressures  from  infin¬ 
ity  to  field  capacity  cover  a  range  of  pore  diameters  up  to  about  8  pm.  Fur¬ 
thermore.  if  one  assumes,  as  did  Campbell  (1988),  that  the  smallest  pores  in 
soil  fill  first  with  water,  then  one  could  always  estimate  the  moisture  content 
below  which  the  soil  structure  can  be  modeled  as  fractal  in  the  following  way. 
First  of  aU,  by  definition 


(  \ 

V 

V 

V 

W 

Y 

T“ 

'J 

=  rt’S 


(15) 


where 

M,  =  volumetric  moisture  content 
=  volume  of  water  in  a  sample 
V,  =  total  volume  of  the  sample 
V,  =  volume  of  voids  in  a  sample 
n  =  sample  porosity 
S  =  degree  of  saturation  of  sample 

But  the  sample  porosity  can  be  rewritten  in  terms  of  sample  dry  density  p^  and 
particle  density  p^,  and  the  degree  of  samration  is  (under  the  pore-filling 
assumptions)  just  the  integral  of  the  pore  distribution  curve  from  zero  volume 
to  the  cumulative  volume  associated  with  the  particular  pore  diameter. 
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The  cutoff  value  of  volumetric  moisture  content  for  fractal  soil  strucmre 
response  can  then  be  written  simply  as 


("J 


'fractal 


(16) 


where  refers  to  the  degree  of  saturation  (or  the  decimal  value  of  cumulative 
pore  volume)  at  the  8-pm  mark  on  a  pore-size  distribution  curve. 

In  summary,  then,  if  one  assumes  that  soil  fills  with  water  beginning  with 
the  smallest  pores,  and  if  one  accepts  the  small-pore  fractal  model,  then  the 
moisture  content  below  which  the  soil  possesses  a  fractal  structure  can  be 
estimated  from  the  sample  dry  density  and  cumulative  pore  volume  given  by 
the  8-pm  pore  diameter  mark  on  the  pore-size  distribution  curve. 


Fractal  Model  Related  to  Particle-Size  Distribution 
Data 

In  the  absence  of  porosimetry  and  pressure  plate  data  (which  is  the  case  for 
this  research),  it  would  be  advantageous  to  estimate  pore  volume  distributions 
from  other  measurements  that  are  more  readily  available. 

At  least  one  simple  model  for  pore  distributions  does  exist;  namely,  the 
physicoemperical  model  offered  by  Arya  and  Paris  (1981).  Their  research 
centered  on  the  need  for  a  model  to  pr^ict  soil  moisture  characteristics  from 
simple  data  like  particle-size  distributions  and  bulk  (dry)  density,  and  their  key 
assumptions  were  as  follows: 

a.  The  particle-size  distribution  curve  can  be  broken  into  segments  cover¬ 
ing  a  small  range  of  diameters,  and  the  solid  fraction  within  each  seg¬ 
ment,  or  particle-size  range,  has  a  bulk  derrsity  equal  to  that  of  the 
natural-structure  sample. 

b.  Each  solid  fraction  consists  of  uniform-size  spheres  having  the  mean 
radius  of  the  fraction. 

c.  The  pores  in  each  segment  are  uniform-size  cylindrical  capillary  tubes 
whose  radii  are  related  to  the  mean  particle  radius  for  that  segment. 

If  there  are  spherical  particles  of  size  in  the  ith  particle-size  segment  of 
the  distribution  curve  (the  following  relationships  assume  that  one  has  a  unit 
mass  of  solid  material:  i.e.,  Z  W-  =  1.0),  then 
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(17) 


I 


where 


=  solid  volume  of  the  ith  fraction  per  unit  sample  mass 

JV,  =  solid  mass/unit  sample  mass  in  the  ith  fraction,  which  comes  from 
the  particle-size  distribution  curve 

fip  =  particle  density 

Arya  and  Paris  then  let  the  total  pore  volume  for  the  ith  fraction  be  repre¬ 
sented  by  a  single  pore  of  radius  and  length  h,  that  threads  the  volume  occ¬ 
upied  by  the  particle-size  fraction.  Thus 


(18) 


where 

=  void  volume  of  the  ith  particle-size  fraction  per  unit  sample  mass 
e  =  void  ratio 

=  volume  voids/volume  solids 

Not  knowing  exactly  how  this  single  pore  threads  the  mass  of  spherical  , 
particles,  Arya  and  Paris  assumed  that  the  total  pore  length  would  be  some¬ 
thing  greater  than  and  introduced  an  empirical  factor  a  defined  by 

h  =  n^lR  (19) 

I  /  / 


Combining  Equations  18  and  19,  one  arrives  at  their  final  relationship 
between  pore  radius  and  particle  radius  in  the  ith  particle-size  fraction 


r.  =  R, 


2en('-“) 


1/2 


(20) 


Chapter  4  Fractal  Geometry  Model  and  Critical  Water  Content 


49 


where  n,  comes  from  Equation  17  and  the  manner  in  which  the  panicle-size 
distribution  curve  gets  subdivided. 

One  now  has  a  model  for  deriving  pore-size  distributions  from  particle-size 
distributions,  densities,  and  an  emperical  parameter  associated  with  soil  texture. 
But  how  can  one  relate  this  new  information  to  soil  moisture  contents?  The 
answer  lies  in  assuming  that  the  smallest  pores  first  fill  with  water  foUowed  by 
larger  and  larger  voids.  Hence,  a  volumetric  water  content  can  be  calculated 
by  summing  the  void  volume  fractions  identified  above: 


9, 


(21) 


where 

Pi  =  the  sample  bulk  density 

At  the  ith  particle-size  fraction,  the  capillary  suction  comes  from  Equation  1 1. 
Thus  one  has  a  model  for  soil  tension  versus  volumetric  water  content  Arya 
and  Paris  exercised  this  model  for  several  soils  and  adjusted  the  empirical 
parameter  a  for  each  to  give  the  best  fit  to  the  experimental  data.  For  five 
different  soils,  the  best  fit  values  of  a  ranged  from  1.35  to  1.39. 

In  summary,  given  some  knowledge  of  the  a  parameter  that  first  appears  in 
Equation  19,  a  model  exists  for  converting  a  particle-size  distribution  curve 
into  a  pore-size  distribution  curve  and  an  opportunity  is  presented  to  apply  the 
fractal  geometry  model  and  compute  a  critical  volumetric  moisture  content, 
below  which  the  soil  structure  can  be  taken  as  a  fractal  geometry. 

As  a  test  of  this  model,  consider  its  reapplication  to  the  data  that  Arya  and 
Paris  provided  in  their  report.  Figures  20  to  24  contain  a  reproduaion  of  their 
particle-size  distribution  curves  for  five  of  the  soils  aixl  soil  mixtures  that  they 
studied  along  with  a  pore-size  distribution  deduced  from  their  soil  tension  data 
using  the  standard  capillary  relationship.  The  contact  angle  was  taken  to  be 
zero  degrees,  and  the  surface  tension  for  pure  water  against  clean  glass  at 
25  °C  was  used  to  convert  pressure  data  to  pore  diameters.  Maximum  volu¬ 
metric  moisture  (100  percent  saturation)  was  taken  to  be  defined  by  the  dry 
density  of  the  sample  as  before.  The  solid  line  on  each  figure  is  the  particle- 
size  distribution  curve,  and  the  dashed  line  is  the  pressure-plate  derived  pore- 
size  distribution.  The  open  circles  on  the  curves  represent  the  results  of 
calculating  a  pore-size  distribution  curve  using  the  model  above. 

Other  than  the  large  pore-size  discrepancy  for  the  sample  shown  in  Fig¬ 
ure  24,  the  two  techniques  for  calculating  pore-size  distributions  (one  from 
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Figure  20.  Comparison  of  measured  and  predicted  pore-size  distributions;  70  percent  silty 
clay,  30  percent  sandy  loam  (measurements  from  Arya  and  Paris  (1981)) 
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Figure  21.  Comparison  of  measured  and  predicted  pore-size  distributions;  loam  40-  to  50-cm 
depth  (measurements  from  Arya  and  Paris  (1981)) 
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Figure  22.  Comparison  of  measured  and  predicted  pore-size  distributions;  40  percent  silty 
clay,  60  percent  sandy  loam  (measurements  from  Arya  and  Paris  (1981)) 
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Figure  23.  Comparison  of  measured  and  predicted  pore-size  distributions;  loam  20-  to  30-cm 
depth  (measurement  from  Arya  and  Paris  (1981)) 
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Figure  24.  Comparison  of  measured  and  predicted  pore-size  distributions;  20  percent  silty 
clay,  80  percent  sandy  loam  (measurements  from  Arya  and  Paris  (1981)) 
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pressure-plate  data  and  one  from  pore-size  distribution  data)  produce  quite 
similar  results.  As  a  brief  aside,  it  is  interesting  to  note  that  Campbell’s 
earlier  assumption  of  median  pore  size  being  equivalent  to  the  median  parti¬ 
cle  size  is  not  true  for  these  data.  One  might  even  say  that  based  on  this 
limited  model  exercise,  the  median  pore  size  is  approximately  an  order  of 
magnitude  smaller  than  the  median  particle  size. 

Now,  how  well  does  the  fractal  model  for  small  pore  structure  apply  to 
the  Arya  and  Paris  data?  From  their  data,  one  can  obtain  the  moisture  con¬ 
tent  for  each  sample  at  field  capacity,  and  from  the  previous  figures,  one  can 
read  5,  ,  the  degree  of  saturation  at  a  pore  diameter  of  8  pm.  Equation  16 
can  then  be  used  to  compute  the  moisture  content  at  the  limits  of  fractal  soil 
structure  (particle  density  is  taken  to  be  2.65  g/cc).  These  results  are 
summarized  in  Table  1  and  show  a  remarkable  correlation. 


Table  1 

Fractal  Cutoff  Moisture  Content  Versus  Measured  Field  Capacity 

Moisture  (data  from  Arya  and  Paris  (1981)) 

Sample 

Dry  Danalty 

s. 

B 

1.400 

0.81 

0.38 

0.39 

C 

1.416 

0.87 

0.41 

0.42 

D 

1.480 

0.67 

0.30 

0.33 

E 

1.456 

0.58 

0.26 

0.28 

F 

1.517 

0.56 

0.24 

0.26 

Fractal  Model  Applied  to  This  Study 

Is  the  fractal  structure  model  compatible  with  the  data  collected  in  this 
study?  Most  importantly,  can  the  electrical  response  of  the  test  soils  be  cor¬ 
related  with  the  fractal  model?  First  of  all,  let  us  state  the  hypothesis; 

a.  Soil  water  fills  the  smallest  pores  first  in  a  dry  sample. 

b.  Soil  structure  can  be  modeled  by  a  fractal  geometry  in  the  range  of 
pore  sizes  in  which  water  is  held  *'v  capillary  forces. 

c.  Within  the  small-pore  fractal  range,  the  fractal  dimension  can  be  cal¬ 
culated  from  knowledge  of  either  the  pore  volume-pressure  relation¬ 
ship  (Equation  14)  or  the  pore  volume-pore  size  relationship 
(Equation  10). 
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d  The  fractal  model  fails  at  pressures  less  than  field  capacity  (1/3  bar), 
which,  using  the  capillary  equation  and  assumed  values  of  surface 
tension  and  contact  angle  for  water,  is  equivalent  to  a  pore  diameter 
of  about  8  (im  (in  other  words,  the  fractal  model  applies  for  pore  dia¬ 
meters  less  than  8  pm). 

e.  The  upper  bound  on  moisture  content  associated  with  the  small-pore 
fractal  model  is  the  field  capacity,  which  can  be  estimated  from  the 
product  of  sample  porosity  (which  is  calculated  from  dry  densitv'  and 
particle  density)  and  the  degree  of  saturation  at  the  8-|im  mark  on  the 
pore  volume-pore  size  curve  (Equation  16). 

/  The  critical  moisture  content  at  which  there  is  a  transition  from  series 
to  parallel-like  electrical  behavior  is  identical  to  the  field  capacity. 


Items  a  -  e  were  examined  in  the  previous  sections.  Item  f  contains  the 
physics  of  the  relationship  between  soil  moisture  levels  and  the  electrical 
response  of  the  soil.  What  it  says  is  that  long-range  connectivity  becomes 
achievable  when  the  pores  governed  by  the  capillary  equation  become  filled 
with  water.  At  moisture  contents  higher  than  this  value,  free  water  exists 
throughout  the  sample;  i.e.,  the  sample  is  no  longer  just  “moist,”  it  is  now 
“wet.” 

The  ideal  set  of  data  required  to  test  this  hypothesis  would  include  enough 
electrical  response  measurements  as  a  function  of  moisture  content  to  deter¬ 
mine  the  critical  moisture  content  and  field  capacity  measurements  taken 
from  pressure  plate  data  at  the  same  sample  densities.  Unfortunately,  these 
data  do  not  exist.  The  best  that  can  be  done  is  to  use  the  electrical  property 
measurements  in  this  study  and  estimate  the  field  capacity  of  the  soil  sam¬ 
ples  from  the  particle-size  distribution  curves  using  the  Arya-Paris  model  to 
generate  a  pore-size  distribution.  This  requires  values  of  the  fitting 
parameter  a. 

Although  Arya  and  Paris  (1981)  were  able  to  use  values  of  the  model 
parameter  that  fell  within  the  range  of  1.35  to  1.39  for  the  five  soils  dis¬ 
cussed  in  the  previous  section,  one  can  argue  that  these  values  may  not  be 
suitable  for  all  soil  samples.  Their  soils  were  well-graded  and  all  contained 
some  sand  as  well  as  some  clay.  The  soils  used  in  this  study  are,  relative  to 
the  Arya-Paris  soils,  poorly  graded.  Recalling  that  the  a  parameter  was  a 
way  of  accounting  for  the  fact  that  the  length  of  a  pore  associated  with  an 
assumed  spherical  particle  is  something  more  than  twice  the  particle  radius,  it 
seems  likely  that  it  could  vary  considerably  for  soils  that  are  poorly  graded 
and  for  which  the  basic  particle  shapes  are  distinctly  different. 

For  example,  pore  spaces  are  likely  to  be  parallel  to  the  flat  plate-like 
structure  of  kaolinite.  The  pore  lengfri  contribution  from  each  particle  should 
be  about  the  particle  diameter.  In  other  words,  a  for  kaolinite  should  be  on 
the  order  of  1 .0.  At  the  other  extreme,  Ottawa  sand  in  a  close-pack 
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geometry  should  possess  pore  channels  that  hug  the  surface  of  the  almost- 
spherical  particles,  which  means  each  particle  contributes  a  path  length  much 
greater  than  its  diameter.  This  would  imply  the  largest  values  of  a. 

Since  there  is  no  a  priori  knowledge  of  the  a  parameter  for  these  soils, 
one  can  take  the  reverse  approach  to  testing  the  above  hypothesis.  Assume 
that  the  critical  moisture  content  is  the  field  capacity,  which,  in  turn,  is  repre¬ 
sented  by  the  degree  of  saturation  shown  on  the  pore  volume  versus  pore 
size  curve  at  the  8-pm  diameter  mark.  Following  the  generation  of  several 
pore-size  distribution  curves  from  the  particle-size  distributions,  estimate  a 
for  each  soil  by  the  value  that  yields  a  pore-size  distribution  that  matches  the 
critical  moisture  and  the  8-pm  intersection.  With  this  value  of  a,  generate 
the  most  likely  pore-size  distribution  curves  for  each  of  the  soils  and  gra¬ 
phically  determine  the  soil  fractal  dimension  for  each  soil  type  using  the 
relationship  in  Equation  10.  The  final  test  of  the  hypothesis  will  be  to  com¬ 
pare  those  values  for  fractal  dimensions  with  data  for  similar  soils  that  are 
published  in  the  open  literature. 

First  of  all.  Figures  25  to  28  show  the  results  of  applying  the  Arya  and 
Paris’  model  to  each  soil  used  in  this  study  with  a  allowed  to  vary  from  1 .0 
to  2.0  in  increments  of  0.1.  Secondly,  values  of  critical  moisture  were 
assigned  to  the  four  soils  based  on  the  moisture-related  data  shown  on  Fig¬ 
ures  11-14.  These  values,  along  with  soil  densities  are  used  to  calculate 
porosities  and  the  required  degrees  of  saturation  shown  in  Table  2.  Armed 
with  this  information,  one  can  return  to  Figures  25  to  28  and  estimate  the 
value  of  a  associated  with  each  soil.  These  values  are  also  shown  in  the 
table. 

Finally,  one  can  now  take  the  slopes  of  the  ca^,'ulated  pore-size  distribu¬ 
tion  curve  for  each  soil  (using  the  estimated  a  value  and  the  fractal  model 
described  above)  and  infer  from  those  results  a  fractal  dimension  for  the  soil 
structure  using  Equation  10.  Figures  29  to  31  show  the  results  of  this  exer¬ 
cise  and  include  the  estimate  of  the  soil  fractal  dimensions,  varying  from  2.0 
for  tan  sand  to  2.9  for  kaolinite.  Ottawa  sand  is  not  reported  because  it 
appears  to  have  no  critical  moisture  electrical  response. 


Supporting  Data 

Although  very  little  data  exist  in  the  literature  on  fractal  dimensions  of 
minerals  and/or  soils,  the  extremes  are  supported  by  at  least  a  couple  of  ref¬ 
erences.  For  example,  Pfeifer  (1984)  calculated  a  fractal  dimension  for  kaol¬ 
inite  of  2.92  using  a  dye  adsorption  and  photometric  measurement  technique. 
On  the  other  hand,  he  and  other  colleagues  (Farin,  Avnir,  and  Pfeifer  1984) 
report  from  another  source  a  dimension  of  about  2.15  for  quartz  sands  using 
a  nitrogen  adsorption  technique. 
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Figure  25.  Predicted  pore  distribution  curves  for  Ottawa  sand 
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Table  2 

Estimating  a  for  the  Soils  in  This  Study 


Soil 

P  (fry 

P  P*rt 

n  j 

H 

Ottawa  Sand 

0.00 

1.7 

2.67 

0.36 

0  00 

1.8* 

Tan  Sand 

0.10 

1.7 

2.66 

0.36 

0.28 

1.8 

Tan  Silt 

0.35 

1.4 

2.71 

0.48 

0.73 

1.3 

Kaolinite 

0.35 

0.9 

2.61 

0.66 

0.53 

0.9 

*  Assumed  to  be  like  that  of  the  tan  sand. 
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Figure  29.  Inferred  fractal  dimension  for  tan  sand 
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Figure  30.  Inferred  fractal  dimension  for  tan  silt 
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Figure  31.  Inferred  fractal  dimension  for  kaolinite 


66 


Chapter  4  Fractal  Geometry  Model  and  Critical  Water  Content 


5  Conclusions 


The  response  of  moist  soils  (or  any  other  moist  heterogeneous  mixture)  to 
active  microwave  sensors  is  controlled  by  surface  geometry  and  the  dielectric 
properties  of  the  medium.  A  review  of  the  literature  reveals  that  what  data 
do  exist  on  dielectric  properties  suffer  from  a  number  of  shortcomings.  In 
most  cases,  the  soils  are  simply  not  properly  characterized  in  terms  of  their 
physical  parameters  such  as  the  distribution  of  particle  sizes,  the  dry  density 
under  which  tests  were  conducted,  or  the  temperature  of  the  soil,  all  of 
which  have  some  impact  on  the  dielectric  properties.  The  volumetric  mois¬ 
ture  content  of  the  soil,  which  is  usually  the  most  important  factor  con¬ 
trolling  dielectric  response,  is  often  not  specified  or  is  not  computable.  And, 
finally,  the  data  are  usually  collected  at  a  limited  number  of  fi^uency  values 
and  over  a  very  limited  range  of  frequencies. 

Attempts  to  model  the  dielectric  response  of  these  soils  are  also  limited  in 
several  respects.  Many  models  are  strictly  empirical  and,  therefore,  apply 
only  to  that  data  set  with  its  few  variations.  Equivalent  circuit  models  have 
been  successfully  applied  to  data  over  a  small  frequency  range,  but  the 
model  parameters  seldom  have  any  physical  meaning. 

This  study  began  with  an  attempt  to  supplement  existing  data  on  the  com¬ 
plex  dielectric  properties  of  moist  soils.  A  coaxial  measurement  apparatus 
for  which  measurements  are  controlled  by  a  vector  network  analyzer  was 
fabricated  to  allow  property  measurements  in  nominally  moist  soils  to  fre¬ 
quencies  above  10  GHz,  which  allowed  one  to  see  the  beginning  of  the 
dipole  loss  mechanism.  Sample  temperatures  were  controlled  by  an  external 
bath,  and  sample  moisture  contents  were  varied  by  incremental  drying  of  the 
samples.  Actual  test  moisture  values  were  determined  after  each  series  of 
measurements  were  completed  by  weighing  the  completely  dry  sample. 
Distilled,  deionized  water  was  used  to  saturate  dry  soil  samples. 

Soil  samples  were  chosen  to  reflect  a  broad  span  of  responses.  They 
varied  from  a  poorly  graded  clean  sand  to  a  very  pure  nonswelling  clay  min¬ 
eral.  The  physical  properties  of  the  soil  samples  were  characterized  by  the 
development  of  grain-size  distribution  curves  and  measurements  of  specific 
gravity  and  specific  surface.  With  hindsight,  either  mercury  porosimetry 
measurements  or  pressure  plate  measurements  should  have  been  conducted  to 
obtain  a  relationship  between  pore  volume  and  pore  size. 
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As  expected,  the  soil  dielectric  response  was  found  to  be  a  strong  function 
of  volumetric  moisture  content  for  nonfrozen  soils.  Low-frequency  losses 
because  of  either  ionic  conductivity  or  Maxwell- Wagner  effects,  or  both, 
were  quite  apparent,  as  were  high-frequency  losses  because  of  the  firee-water 
dipole  relaxation  loss  mechanism.  Minimum  losses  occur  in  the  1-  to  2-GHz 
range. 

Sample  temperature  did  not  have  a  significant  effect  on  the  dielectric 
response  of  nonfrozen  soils  over  a  range  of  temperatures  from  0  to  40  °C  at 
any  frequency.  Data  collected  on  frozen  soils  at  -10  and  -5  °C  also  did  not 
show  any  significant  difference.  The  pronounced  transition  from  frozen  to 
nonfrozen  soils  took  place  within  the  -5  to  0  °C  range,  but  the  spacing  of 
data  within  this  range  was  not  enough  to  completely  define  that  transition. 

Data  collected  in  the  study  supported  the  claim  of  previous  researchers 
that  there  is  a  definite  transition  in  normal  soils  from  a  series-like  capacitive 
response  to  a  parallel-like  response  with  increasing  volumetric  moisture  con¬ 
tent.  For  the  silt  and  clay  samples,  that  transition  occurred  in  the  30-  to 
3S-percent  moisture  content  range,  while  for  the  well-graded  sand,  it 
occurred  in  the  10-  to  15-percent  range.  For  equal  moisture  contents,  then, 
the  sands  had  a  higher  permittivity  than  silts  and  nonswelling  clays. 

Two  approaches  were  taken  to  modeling  the  dielectric  response  of  these 
soils.  One  was  an  equivalent  circuit  that  allowed  for  both  series-like  and 
parallel-like  responses  of  soil  and  water  elements.  An  iterative  solution  tech¬ 
nique  was  used  to  arrive  at  the  weighting  parameters  for  each  element  that 
varied  in  a  sensible  way  with  moisture  content.  The  permittivities  assigned 
to  the  circuit  elements  were  held  fixed  at  acceptable  values,  and  excellent  fits 
to  the  data  over  the  entire  range  of  frequencies  were  achieved  for  ail  of  the 
soil  types. 

The  second  model  was  one  that  assumed  the  small  pore  geometry  of  a 
porous  medium  like  soils  could  be  described  by  fractal  methods.  The  result 
was  an  expression  for  pore  volume  as  a  function  of  pore  size  or  capillary 
pressure  that  included  a  fractal  dimension  of  the  soil  structure.  When 
applied  to  a  set  of  pressure  plate  data  taken  from  the  open  literature,  this 
model  showed  that  the  fractal  geometry  assumption  held  for  the  small-pore 
portion  of  the  data  and  collapsed  at  the  field  capacity  of  the  soils  (1/3  bar 
pressure,  or  8-pm  pore  diameter).  The  model  goes  on  to  hypothesize  that  the 
critical  moisture  content  of  the  soils  at  which  its  electrical  response  makes  a 
transition  from  series  to  parallel  behavior  is,  in  fact,  the  moisture  content  at 
field  capacity.  Lacking  a  measured  pore  volume  distribution  for  these  soils 
in  terms  of  either  pressures  or  pore  size,  a  pore  volume-pore  size  distribution 
was  estimated  from  the  measured  particle  size  distribution  curves  and  used  to 
evaluate  bounds  on  the  small-pore  fractal  dimension  of  the  idealized  soils 
used  in  this  study.  The  results  indicate  that  while  real  soils  are  likely  to  have 
fractal  dimensions  of  2.5  to  2.7,  the  value  for  coarse  sands  tends  toward  2.0 
and  the  value  for  fine  nonswelling  clay  minerals  tends  toward  3.0. 
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Appendix  A 
Soil  Properties 


Particle-size  distribution  curves  for  titch  of  the  four  soils  used  in  t.iese 
experiments  were  developed  by  technicians  i"  the  Geotechnical  Laboratory  at 
the  U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS,  and 
are  shown  on  the  following  four  figures.  Also  included  on  each  figure  is  the 
specific  surface  for  that  soil  that  was  measured  by  laboratory  personnel  at 
Quantachrome  Corporation,  in  Syosset,  NY,  using  a  nitrogen  gas  adsorption 
technique  (Krypton  gas  for  the  Ottowa  sand). 
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Gradation  curve  for  Ottawa  sand 


Figure  A2.  Gradation  curve  for  tan  sand 
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Figure  A3.  Gradation  curve  for  tan  silt 


Appendix  A  Soil  Properties 


A5 


Figure  A4.  Gradation  curve  for  kaolinite 


Appendix  B 

Fractal  Models  of  Soil  Structure 


Fractals 

In  recent  years,  some  very  interesting  woric  has  been  done  on  the  conduc¬ 
tive  behavior  of  porous  media  using  the  concept  of  self-similar  geometries, 
or  fractals,  to  model  the  behavior  of  certain  physical  quantities  of  that  media. 
Mandelbrot,  the  guru  of  fractals,  defines  self-similar  shapes  as  those  in  which 
a  certain  part  of  the  shape  can  be  broken  up  into  N  smaller  parts,  each  look¬ 
ing  like  the  original,  but  reduced  in  size  by  a  fractional  factor  labeled  r 
(Mandelbrot  1983).'  The  fractal  dimension  D  of  this  self-similar  shape  can 
be  written  as 


D  = 


log(A^) 

•ogf- 


(Bl) 


or 


N  =  I 


(B2) 


Mandelbrot  insisted  that  the  exponent  D  be  thought  of  as  a  dimension, 
because  it  arose  from  a  method  of  measuring  the  perimeter  of  an  object  hav¬ 
ing  an  irregular  boundary. 

Fractals  are  often  used  to  generate  complex,  or  textured,  curves  or 
surfaces  from  very  simple  initiator  geometries.  For  example,  consider 
Figure  1  in  which  the  initiator  (or  the  initial  segment  of  the  geometry  that 


'  Refemices  cited  in  this  appendix  are  located  at  the  end  of  the  main  text 
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Figure  B1.  A  fractal  snowflake  of  dimension  1.5  (from  Mandelbrot  (1983)) 


will  be  fractalized)  is  a  solid  square  (precisely  stated,  each  side  of  the  square 
is  an  initiator),  and  the  generator  (or  the  desired  geometry)  is  the  broken-line 
segment  as  shown.  Moving  clockwise  from  the  initiator,  the  first  two  appli¬ 
cations  of  the  generator  are  shown  at  full  scale,  while  the  next  two  are 
enlarged  to  show  the  edge  detail.  Figure  B2  shows  another  construction 
initiated  with  a  solid  square  in  which  the  generator  results  in  the  formation  of 
both  “lakes”  and  “islands.” 


Soil  Structure 

A  third  example  of  a  fractal  construction  that  is  useful  for  studying 
porous  media,  and  may  be  particularly  useful  for  studying  soil  structure,  is 
that  shown  in  Figure  B3  in  which  the  initiator  is  a  cube  with  square  holes 
centered  on  each  face  and  joining  in  die  center.  The  generator  is  the  same 
object  but  reduced  by  a  factor  of  1/3.  Therefore,  in  the  first  stage,  the 
initiator  is  divided  into  20  smaller  cubes,  resulting  in  a  fractal  dimension  of 
2.7268.  Mandelbrot  referred  to  this  construction  as  a  Menger  sponge. 

Friesen  and  Mikula  (1987)  used  the  Menger  sponge  to  model  the  fractal 
dimension  of  the  surface  area  of  coal  in  such  a  way  that  the  slope  of  the 
isotherm  (fractional  volume  change  versus  pressure  for  porosimetry  measure¬ 
ments)  plotted  on  a  log-log  scale  is  linearly  related  to  the  fractal  dimension. 
At  least  one  recent  study  attempted  to  measure  the  fractal  dimensions  of 
sandstone  solid-void  interfaces  (Katz  and  Thompson  1985).  Measured  values 
of  interface  fractal  dimensions  on  five  different  samples  using  electron 
microscope  techniques  varied  from  2.57  to  2.87. 
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Figure  B2.  Fractal  islands  and  lakes  of  dimension  1.6131  (from  Mandelbrot  (1983)) 

Work  has  been  done  very  recently  on  modeling  the  microstructure  of  soil 
using  fractals  (Krepfl,  Moore,  and  Lee  1989;  Moore  and  Krepfl  1991).  A 
two-dimensional  representation  of  a  three-dimensional  model  using 
hexagonal-shaped  flakes  as  the  generator  is  shown  in  Figure  B4.  Such  mod¬ 
els  may  be  useful  as  a  representation  of  pure  clays.  This  particular  fractal 
structure  is  strongly  reminiscent  of  microphotographs  of  the  pure  kaolinite 
mineral  structure. 
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